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ABSTRACT 

We present new high-resolution near-IR spectroscopy and OH maser obser- 
vations to investigate the population of cool luminous stars of the young massive 
Galactic cluster RSGCl. Using the 2.293/im CO-bandhead feature, we make 
high-precision radial velocity measurements of 16 of the 17 candidate Red Super- 
giants (RSGs) identified by Figer et al. We show that F16 and F17 are foreground 
stars, while we confirm that the rest are indeed physically-associated RSGs. We 
determine that Star F15, also associated with the cluster, is a Yellow Hypergiant 
based on its luminosity and spectroscopic similarity to p Gas. Using the clus- 
ter's radial velocity, we have derived the kinematic distance to the cluster and 
revisited the stars' temperatures and luminosities. We find a larger spread of 
luminosities than in the discovery paper, consistent with a cluster age 30% older 
than previously thought (12±2Myr), and a total initial mass of (3± 1) x lO^M©. 
The spatial coincidence of the OH maser with F13, combined with similar radial 
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velocities, is compelling evidence that the two are related. Combining our re- 
sults with recent SiO and H2O maser observations, we find that those stars with 
maser emission are the most luminous in the cluster. From this we suggest that 
the maser-active phase is associated with the end of the RSG stage, when the 
luminosity- mass ratios are at their highest. 

Subject headings: open clusters & associations, supergiants, stars:evolution, starsdate- 
type, masers 



Introduction 



The Red Supergiants (RSGs) repres ent a key evolutionary p hase in the life-cycle of 
stars with initial masses of ~8-3OM0 (e.g. iMeynet fc Maederll200Cll ). Though comparatively 
brief, the mass-loss rate in this stage can be many orders of ma gnitude greater than on the 
main-sequence (e.g. iRepolust et al.l l2004l : Ivan Loon et al.ll2005l ). and the mass lost in the 



RSG phase can determine the terminal mass of the st ar, the appearanc e of the supernova 
(SN) explosion, and nature of the stellar remnant (e.g. iHeger et al.l 120031 ). 



The study of RSGs is hampered by low-number statistics - until recently only ~200 
were known in the Galaxy, and around 40 in the Large Magellanic Cloud (LMC). Further, 
the majority of these stars are isolated, and hence have a variety of ages, initial masses and 
metallicities. Ideally, we would like to study large numbers of RSGs coeval clusters, where 
we can be confident that the variables of metallicity and initial stellar mass are fixed. 

Two recent discoveries now present us wit h the oppor t unity to study statistically- 
significant numbers of RSGs in such clusters. In iFiger et al.l (l2006l. hereafter FMR0 6) the 
discovery of an apparent cluster of 14 RSGs was presented; while IPavies et aL ( 2007 . here- 
after DFK07) report on the remarkable stellar population of a second cluster in the same 
region, which it was shown contains 26 RSGs. In order to use these two clusters as testbeds 
with which to study the pre-SN evolution of massive stars, we must first quantitatively study 
the properties of the stars within each. 

In the case of the second cluster, hereafter RSGC2 (also known as Stephenson 2), 
DFK07 used high-resolution spectroscopy to obtain accurate radial velocities of many stars 
in the region of the cluster, and were able to separate proper cluster-members from back- 
ground/foreground stars. Further, from the radial velocities they made a quantitative dis- 
cussion of the kinematic distance to the cluster, thus enabling them to determine the stars' 
luminosities, temperatures, ages and initial masses. 
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The details of the stellar properties of the first cluster, hereafter RSGCl, are less well 
constrained. In the discovery paper FMR06 showed with low-resolution spectra that the stars 
were all late-type; but with no radial velocity information, they relied on the stars' similar 
near-IR colours to argue that the stars were all at the sa me distance. They deter mined this 



distance by associating a nearby OH maser, detected by lBlommaert et al.l (119941 ). with the 
RSGs. However this maser was sin gly- peaked - such masers are formed in the winds of RSGs 
and hence are typically doubly-peaked, with separation twice the outflow speed. If the second 
peak was missed in the original observations, this would lead to an incorrect radial velocity 
for the cluster, and hence kinematic distance. Consequently, the cluster membership, as well 
as the properties of the luminous cool stars of RSGCl are poorly-constrained when compared 
to RSGC2. 

Here, we study the cool supergiants of RSGCl with data of similar quality to that pre- 
sented in DFK07 for RSGC2. We present high resolution spectroscopy of the CO bandhead 
feature of 16 of the 17 i^-bright stars in the field identified in FMR06, allowing us to deter- 
mine accurate radial velocities of the stars, and establish the cluster membership. We also 
present new observations of the OH maser source OH25.25-0.16, showing that it is indeed 
doubly-peaked, and that the central radial velocity of the profile is consistent with the av- 
erage stellar radial velocity. We use this value, in conjunction with contemporary Galactic 
rotation curve parameters, to reappraise the distance to the cluster, the stars' temperatures 
and luminosities, and the age and intial mass of the cluster. 

We begin in Sect. [2] with a description of the observations and data reduction steps, 
and we describe the results and analysis of the data in Sect. [31 In Sect. S] we derive the 
cluster's age and mass, and compare with similar analyses of RSGC2. Finally, we use the 
stellar population of RSGCl to investigate the maser-active phase of RSGs. 



2. Observations & data reduction 

2.1. Radio observations & data reduction 

OH maser observations at 18 cm were carried out with the VLA0 on 25 May, 2006 in 
the AB configuration. The 256-channel spectrum had a bandwidth of 1.56 MHz giving a 
velocity resolution of 1.1 km s~^ and was centered at wlsr = 100 km s~^. The observation 
consisted of a single pointing with J2000 coordinates of (18h37m52s,-06°53'40") and had a 



^The Very Large Array (VLA) is operated by the National Radio Astronomy Observatory under cooper- 
ative agreement with the National Science Foundation. 
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Table 1:: Data for the stars observed. Designations are from FMR06, coordinates and apparent 
magnitudes are from 2MASS. 



ID 




a 

J2000 


5 






niKs 


FOl 


18 


Q7 
O ( 


00. zy 


p. 

-0 


52 32 2 


9.748 


6.587 


4 




F02 


18 


37 


55.28 


P. 

-0 


52 48.4 


9.904 


6.695 





uzy 


F03 


18 


37 


59.72 


-6 


53 49.4 


9.954 


6.921 





QQQ 

ooo 


F04 


18 


37 


50.90 


-6 


53 38 2 


9.658 


6.803 







F05 


18 


37 


55.50 


-6 


52 12.2 


10.547 


7.178 





Ooo 


F06 


18 


37 


57.45 


-6 


53 25.3 


9.866 


7.038 




613 


F07 


18 


37 


54.31 


-6 


52 34.7 


9.941 


7.065 


5 


631 


F08 


18 


37 


55.19 


-6 


52 10.7 


10.772 


7.330 


5 


654 


F09 


18 


37 


57.77 


-6 


52 22.2 


10.262 


7.240 


5 


670 


FIO 


18 


37 


59.53 


-6 


53 31.9 


10.179 


7.218 


5 


709 


Fll 


18 


37 


51.72 


-6 


51 49.9 


10.467 


7.325 


5 


722 


F12 


IJ 


, 38 


3.30 


-6 


52 45.1 


10.143 


7.238 


5 


864 


F13 


18 


37 58.90 


-6 


52 32.1 


10.907 


7.716 


5 


957 


F14 


18 


37 47.63 


-6 53 2.3 


10.495 


7.576 


6 


167 


F15 


18 


37 57.78 


-6 52 32.0 


10.651 


8.070 


6 


682 


F16 


IJ 


, 38 


1.30 


-6 52 52.0 


13.617 


9.608 


7 


558 


F17 


18 


37 48.77 


-6 53 7.7 


12.763 


10.188 


9 


003 



field of view of about 27'. The beam size is 3.5" x 1.8". The total integration time was 340 
minutes, giving an rms error of about 4 mJy beam~^ per 4 kHz channel. 

The data were calibrated using the AIPS package. The flux and bandpass calibrator was 
1331+305, while 1822-096 was used as a phase calibrator. The phase calibrator was mildly 
resolved in this observation, so baselines with uv distances greater than 50 kX were removed 
when solving for a phase solution. Minimal flagging was required to remove interference. 
The astrometric accuracy of our positions is dominated by the uncertainty in the position of 
our phase calibration source, which is 0.3". The compact source GPSR5 25.266-0.161 has a 
J2000 position of (18h37m57.9934s, - 06°53'30.973"), con sistent with previous measurements 



with an astrometric accuracy of 2" ( iBecker et al.lll994l ). 
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2.2. High-resolution spectroscopy 

2.2.1. Observations 

Observations were taken with NIRSPEC, the cross-dispersed echelle spectrograph mounted 
on Keck-II, during the night of S*'^ May 2006. We observed stars FOl - F17, with the ex- 
ception of F12 which was missed due to time constraints. F15 was re-observed in a separate 
observing run on 12*'^ Aug 2006, using the same technical setup as described below. Table 
[T] lists the coordinates and 2MASS magnitudes of the 17 stars, as identified in FMR06. The 
stars' locations are illustrated in Fig. [2l 

We used the NlRSPEC-7 fiher and 0.576"x24" slit. Setting the dispersion angle to 
62.53°, and cross- disperser angle to 35.53°, this gave us a spectral resolution of ~17,000 in 
the wavelength range 1.9-2.4/im. We integrated for 20s in two nodded positions along the 
slit for each star. We observed the BOV star HD 171305 as a telluric standard. We used a 
continuum lamp to obtain fiat-field frames, while for wavelength calibration, we observed Ar, 
Ne, Xe and Kr lamps to obtain as many spectral lines as possible in the narrow wavelength 
range. To sample the gaps between these lines, we also observed the continuum lamp through 
an etalon filter. 



2.2.2. Data reduction 

To remove sky emission, dark current and bias level, we subtracted nod pairs of spec- 
tra. Fluctuations in pixel-to-pixel sensitivity were corrected for by dividing through by the 
normalized fiat-field frame. 

The optics of NIRSPEC produce spectral orders which are warped in both the spatial 
and dispersion directions, and before the spectra can be extracted this warping must be 
corrected for in a process known as rectification. The methodology we use is the same as 



that in DFK07, and is explained in more detail in iFiger et al.l (120031 ). Here we give a brief 
summary. 

Spatial rectification is done by adding nod-pair spectra and fitting the spectral traces 
with a polynomial. Rectification in the dispersion direction is more complicated; it involves 
obtaining accurate wavelengths of the etalon lines using the arc frames, and assuming that 
the wavelength of the nth order etalon-line A„ is governed by the separation of the etalon 
plates t via the relation A„ = t/2n (the etalon equation). 



The arc lines were used to get initial estimates of the etalon-line wavelengths, and hence 
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of the etalon plate separation. The etalon-hne wavelengths were then recomputed using 
the etalon equation, and used to re-estimate the wavelengths of the arc lines. The etalon 
thickness was fine-tuned in an iterative process until the residuals between the measured and 
predicted arc-line wavelengths across all orders were minimized. 

After rectification, the spectra were extraced from each frame by summing the pixels 
across the trace in each channel. Shifts between spectra of up to 4km s~^ (;^1 pixel), caused 
by the star not being quite in the centre of the slit, were corrected for by cross-correlating 
the atmospheric CO2 feature at 2.05/im in each spectrum. 

The accuracy of the final wavelength solution is determined from the residuals between 
the observed and predicted arc-line wavelengths in the etalon-fitting process described above, 
and is better than ±4 km s~^. The internal error between spectra, from the CO2 telluric 
feature, is <Clkm s~^, and so is dominated by systematics in our analysis process which we 
estimate to be ±lkm s"^ (DFK07). 

3. Results &: analysis 
3.1. The maser source OH25. 25-0.16 



As mentioned in Sect. [T], the 1612MHz OH maser forms in the outflows of RSGs, far 
above the stellar surface. The velocity profiles are therefore typically doubly-peaked, with a 
separation twice the termin al velocity of the outflow centred on the star's systemic velocity. 
However, when observed by lBlommaert et al.l (119941 ). OH25.25-0.16 appeared only as a single 
peak with a radial velocity wlsr = 102.2km s~^. 

Our new observation of the OH maser source is shown in velocity-space in Fig. [H Here 
it can be seen that we clearly detect the second peak. The fiux-weighted mean velocities of 
the two peaks are 103.8 ± 0.1km s~^ and 138.1 ± 0.9km s~^, calculated using all channels 
with emission greater than Sa above the background. The average velocity of the peaks is 
120.9 ± 0.9km s"^, consistent w ith the average velociti e s of t he SiO masers in the cluster, 
120.7 ± 3.2 km s~^, observed by iNakashima &: Deguchi J2006h. The implied o utflow speed 
is 17.1±0.6km s~^, a typical outflow speed for RSGs (Richards &: Yateslll998l ). and similar 



to the outflow speed of S Per (~16km s ^ Diamond et al. 79871). which occ upies a similar 



location in the HR-diagram as the stars of RSGCl (iGahm fc Hultqvistlll976l . see Sect. Hj). 

In Fig. [2] we overlay a contour plot of the 1612MHz emission on the 2MASS J^^-band 
image of the cluster. We find the positional centroid of the OH maser to be 18'*37™58.882*, 
-06°52'32.28" (J2000), with a positional uncertainty of 0.3". The J2000 position of the 
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Fig. 1. — : Spectrum at 1612MHz of the source OH25.25-0.16. Our new observations clearly show 
the second peak at ~135km s^^. Errorbars show ztlo" uncertainty and triangles show Scr upper 
limits in the flux density. 

maser is consistent with previous measurement by Blommaert et al. (1994), which had a 
positional accuracy of 4". The maser is also spatially coincident with star F13, whose 
2MASS coordinates are 18'^37™58.908^ -06°52'32.11" (J2000), with positional uncertainty 
0.06". 

The c entral radial velocity of the O H maser is also consistent with the SiO observation 
of F13 bv iNakashima fc Deguch] j20od lR. vlsr,fi3 = 120.5 ± 2.0 km s-\ and with the CO- 
bandhead radial velocity measurement of F13 presented in this paper, 125.4±4km s~^ (see 
Sect. 13.2.11) . From this evidence, it seems highly likely that the OH maser originates in the 
outflow of F13. 

In addition to the OH maser, we also detect the continuum sources GPSR5 (25.266- 
0.161, 25.25 2-0.139, 25.237-0.15 ) . The se sources were shown conclusively to be extragalactic 
in origin by iTrejo fc Rodrfguea (120061 ) . 



^This radial velocity was found from the average velocity of the high- and low- velocity edges. The centroid 
of the peak of this source was found to have a radial velocity of 116.5±2km s^^ 
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Table 2:: Observed data for RSGs. (1): the stellar IDs (fro m FMR06); (2): the radial ve locity of 
each star (ib4km s~^), and the radial velocity measured by iNakashima &: Deguchil ( 20061 ) . where 
available (±2 km s~^); (3): effective temperature; (4): spectral type, accurate to ±2 subtypes; 
(5) derived extinction towards each star; (6): absolute K-hand magnitude; and (7): bolometric 
luminosity. 



(1) 

ID 


(2) 

lisR (ND06)i 
( km s~^) 


(3) 

Teff (K) 


(4) 

Spec Type 


(5) 




(6) 
Mk 


(7) 


FOl 


129.5 


(117.7) 


3450±127 


M5 


2.58±0 


09 


'^-0.30 


^•^^-0.13 


F02 


114.2 


(119.7) 


3660±127 


M2 


2.83±0 


07 




'-'•^"-0.13 


F03 


127.2 




3450±127 


M5 


2.46±0 


09 


-11.28;°J^ 




F04 
F05 
F06 


121.2 
124.8 
120.7 


(124.3) 


3752±117 
3535±130 
3450±127 


Ml 
M4 
M5 


2.46±0 
2.77±0 
2.19±0 


04 
08 
09 


-11 24+°-32 
-11 S6+°-3^ 
-10.70t°i^ 


cr qo+O.ll 
'-'■'-'^-0.13 

S 29+0-12 
"-••^^-o.u 

5.oo;°:}^ 


F07 


121.6 




3605±151 


M3 


2.33±0 


12 


-io.8it°:i 


5.10t°il 


F08 


128.2 




3605±151 


M3 


2.84±0 


12 


1 1 qq+0.36 
-ii.OO_Q 32 


5 30+0-13 


F09 


121.6 




3399±150 


M6 


2.44±0 


08 


-10.92^°:i 


c- Q7+0.12 


FIO 


122.0 




3605±151 


M3 


2.45±0 


12 


-10.86^°^^ 


c; 10+0.13 
^•-■-^-0.14 


Fll 


124.1 




3535±130 


M4 


2.63±0 


08 


-1 1 ns+°-34 

ii-UO_Q 3Q 


5 16+0-12 


F13 
F14 


125.4 
122.0 


(120.5) 


4015±140 
3605±151 


K2 

M3 


3.19±0 
2.29±0 


09 
12 


-11 39+0-34 
-10.25+°:^^ 


5.45;°:- 

4.87;°:}! 


F15 


120.8 




6850±350 


GO 


2.65±0 


04 


1Q Q7+0.40 


5.36^°:}^ 


F16 


42.6^ 
















F17 


33.2^ 

















^We quote the Nakashima fc Deguchi ( 2006h values measured by taking the average of the high- and low- 
velocity edges of the maser profiles. 

^Stars F16 & FIT are determined to be foreground stars, and due to the uncertainties in reddening and 
distance we derive no stellar parameters for these stars. 
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Fig. 2. — : Contour-plot of the 1612Hz radio image, overlayed on the 2MASS K^-band image of 
the same region. The image is centered on the indicated coordinates (epoch J2000). Contours are 
drawn at 5, 10, 20 and 40(T above the background; the \a level corresponds to 0.87mJy/beam. 
The position of the OH maser is indicated, and is coincident wit h star F13. The other ide ntified 
radio sources in the field, labelled in the figure, were shown by iTreio &: Rodriguea (j2006l ) to be 
extragalactic in origin. The stars FOl to F17, as designated by FMR06, are also labelled. 



3.2. The high-resolution spectra 

The high-resolution observations of the region around the CO-bandhead feature at 
2.293/im for the 17 K-bright stars are shown in Fig. [3l All stars observed, with the ex- 
ception of F15, show the feature strongly in absorption; F15 has it weakly in emission. 
Quantitative analysis of the high-resolution spectroscopy results are described below. 



- 10 - 




2.29 2.30 2.31 

Wavelength (wm) 



Fig. 3. — : High-resolution spectra of all stars observed in the region of the CO-bandhead feature. 
The dotted blue line indicates the zero velocity of the blue edge of the CO bandhead. 
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3.2.1. Stellar radial velocities 



The radial velocities of each star serve two purposes; firstly they allow us to distinguish 
between genuine cluster stars and foreground stars with similar colours; and secondly they 
allow us to derive a kinematic distance to the cluster. 

The differences in radial velocities of the stars can be seen qualitatively in Fig. [3] - the 
blue-edges of the CO bandhead in stars F16 and F17 are noticably blue-shifted compared to 
stars FOl - F15, which all have very similar radial velocities. 

In order to accurately quan tify the rad i al velo cities of the stars, we implemented the 
same technique as presented in iFiger et al.l (120031 ) and DFK07. We cross-correlated the 
spectra shown in Fig. [3| w ith that of the high-resolution spectrum of Arcturus presented in 
Wallace fc Hinkld (Il996a( ). which had been degraded to the same spectral resolution as our 
data. For star F15, which has CO in emission, we inverted the spectrum before analysis. 
We experimented with isolating different spectral ranges during the cross-correlation, to 
determine the robustness of our measurements. We found that the measured velocities 
were stable to within ±1 km s^^, which we take to be the internal error between individual 
measurements. The absolute uncertainty on the measurements is limited by the accuracy of 
the wavelength solution, ±4 km s~^ (see Sect. [2]) • 

We find that stars FOl - F15 all have radial velocities in the range wlsr ~1 15-125 km s~^, 
while stars F16 and F17 have wlsr of 33 km s~^ and 43 km s~^ respectively. From this, we 
conclude that the two faintest stars observed are foreground stars, while stars FOl - F15 
are physical members of the cluster. As F12 was unobserved, the status of this star is still 
unclear. We note that the CO absorption strengths of stars F16 and F17 are much lower than 
those of the RSGs, and instead are more typical of less luminous stars. This is consistent 
with these stars being foreground objects. 

The measured radial velocities of all stars observed are listed in Table [H Also listed in 
Table [H where availab l e, are the radial velocities determined from SiO maser emission by 
Nakashima fc Deguchil (120061 ). We see that the two measurements of stars F02, F04 and F13 
are within ~2(j, while for FOl it is ^3cr. A s SiO masers are commonly thought to trace the 



stellar systemic velocity (jjewell et al.lll99ll ). we would ex pect the two velocity measurem ents 
to agree well. However, we note that the observations of iNakashima fc Deguchil (120061 ) may 
have been hampered by their beam-size; their observations of FOl, F02 and F13 have at least 
one other RSG within the beam FWHM. In the case of F04, which is well separated from 
the other RSGs, the CO bandhead and SiO maser measurements are in excellent agreement. 
As mentioned in Sect. 13.11 the measured radial velocity of F13 is in excellent agreement with 
observations of the OH and SiO maser sources at the same location. 
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3.2.2. Distance to RSGCl 



We take the mean radial velocity of the stars observed at high spectral resolution and 
compare with the Galactic r otation curve, using the contemporary measurements collated by 
Kothes fc Doughertyl (120071 ). The mean radial velocity of the stars from the CO observations 
is 123. Oil. km s~^, with an uncertainty determined from Poisson statistics of the measure- 
ments. The absolute uncertainty on the radial velocity is therefore dominated by that in the 
wavelength solution, ±4 km s ~^. This compares well to the average radial velocity found by 
Nakashima k Deguchil (120061 1 . ~120±2km s'^ from their SiO maser observation^, and the 
central velocity of our new 1612MHz OH maser observation, 120.9 ± 0.9 km s^^. 

In Fig. m we compare our radial velocity to the Galactic rotation curve in the direc- 
tion of / = 25.15°, b = —0.15°. W e use the distance to the Galactic centre -Dgal = 
7.5 =fc 0.3kpc (lEisenhauer et al.l 120051). and solar r otational velocity Qq = 214 ± 7 km s~^ 



(IFeast fc Whitelock 



1997 



Reid &: Brunthalerl 120041 ) . We use the uncertainties on these val- 



ues to construct 'maximal' and 'minimal' rotation curves in Fig. HI 

The cluster radial velocity actually extends beyond the asymptotic point of the curve; 
however it lies well within the two 'error' curves, and therefore could simply be due to 
the uncertainties in -Dgal and 9©. To determine the distance to the cluster, we take the 
average of the two points where the radial velocity intercepts the 'maximal' rotation curve. 
This gives us a kinematic distance to the cluster of 6.60±0.89 kpc, slightly larger than the 
distance quoted in FMR06 when using the radial velocity of the singly-peaked OH maser. 



3.2.3. Effective temperatures 



To determine the spectral-types of the RSGs, and hence t h eir ef fective temperatures, 
we used the same empirical method described in iDavies et al.l (120071 ). We compared the 
equivalent width (IVa) of the CO bandhead absorption with that of tem plat e stars, taken from 
the ca talogues of lKleinmann fc Halll (119861 ). IWallace k. Hinkld (Il996bl ) and I Wallace fc Hinkle 
( 119971 ). We defined a measurement region of 2.294-2.304/im, and defined the continuum as 
the median of the range 2.288-2.293/im. We estimated the uncertainty by repeating the 
measurements with slightly adjusted continuum regions, and found that measurements were 
stable to ~1A, or ~5%. To find the spectral-types of the RSGs, we compared the W\ 
measurements with a linear fit to the W\ of the template stars as a function of spectral- 



•^These values are the average of their two different methods of measuring the radial velocity of each star 
from the maser line profile 
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Fig. 4. — : Galactic rotati on curve in the direct i on of RSGCl ([Brand &: Blitzj Il993l ). using the 
measurements of collated by iKothes h Doughertyl (j2007l ) . 



type. Using this method, we are able to determine spectral-types to within ±2 subtypes 
(see DF K07). In conver t ing sp ectral-type to effective temperature, we used the temperature 



scale of iLevesque et al.l (120051 ) . The derived spectral- types and effective temperatures for 



the RSGs are listed in Table El 

For star F15 the method described above breaks down, as this star has CO in emission. 
The star's radial velocity suggests that it is part of the cluster, and hence a supergiant. It was 
assigned the spectral-type G6 I in FMR06, and hence deemed to be a 'Yellow Hypergiant' 
(YHG), based on its weak CO-bandhead absorption. Given that this was only a marginal 
detection of CO absorption here we reappraise the star's spectral type. 

Figure [5] shows the spectra of F15 on the two occasions it was observed. It can be seen 
that the CO-bandhead, which was present in emission in May, was not detected at all in 
August. Over the rest of the star's spectrum, no discernable variability is observed. 
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For comparison, Fig. [5] shows similar template spectra of 7 Cyg, spectral-type F8 lab, 
and a Sge, spectral-typ e Gl III, which have be en resampled to the same spectral resolution 
as our data (taken from [Wallace fc Hinkldl 19971 ). Though admittedly a Sge has a lower lumi- 
nosity class than that inferred for F15, it serves to give some insight into F15's temperature. 
At spectral-type Gl, the CO-bandhead absorption is still seen, albeit weakly. The star also 
shows the atomic absorption lines of Al I, Fe I, Mg I, and Br 7. However, at this temperature 
molecular absorption, mostly from CN, can be seen as the undulating 'noise'-like features 
throughout the spectrum. At spectral-type F8, the CN and CO absorption are gone, while 
the atomic absorption lines remain. 

From these comparison spectra, we assign a spectral-type to F15 of GO, ±2 subtypes, 
and hence an effective temperature of Tgfj = 6500 — 7200K. Indeed, the star is very similar 
spectroscopically to p Gas, one of the archety pal YHGs, which a lso shows transient GO 
emission / absorption (see spectra presented in lGorlova et al.ll2006l ). 



3.2.4- Extinction 

The extinction towards each star is measured by comparing the 2MASS infrared colours 
of the stars to the intrinsic colours of supergiants with the same spectral-typ e. For the 



RSGs we use the observations of Galactic stars wit h luminosity clas s lab from lElias et al. 



(119851 ): while for F15 (the YHG) we use the tables of lKoornneeil ( 1l983l ). We define the excess 



between the observed and intrinsic colours as the reddening towards eac h star. We convert 



this re ddening to an extinction towards each star using the relationship of iRieke fc Lebofsky 



(]1985f ) 



A 



E 



X-Ks 



where A is the wavelength appropriate for the 2MASS J or H filters. To determine the 
uncertainty in each extinction measurement we derive extinctions for the upper and lower 
limits to each star's spectral-type. Where the difference between the derived A^giJ — Ks) 
and Aks{H — Ks) values is outside this uncertainty, we adopt half this difference as the error 
in the measurement. 

Using this method, we find a median extinction of Aks{J — Ks) = 2.58, and AxgiH — 
Ks) = 2.62. Each have uncertainties of 0.07mags from Poisson statistics, and are therefore 
in good agreement with one another. We adopt the mean of these measurements, A^s = 
2.60±0.07 as the median extinction towards the cluster. We note that the extinction towards 
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the YHG, v4k,fi5 = 2.65 ± 0.02, is consistent with that derived for the RSGs. 

The median cluster extinction is shghtly lower than the A^g = 2.74 ± 0.02 derived 
in FMR06. We consider the latest measurement to be the more reliable, due to the extra 
assumptions used in FMR06: instead of dereddening each star according to the intrinsic 
colours appropriate for its spec t ral- ty pe, they dereddened all stars to the mean colour of 



M supergiants from lElias et al.l (119851 ). The ranges in colours of M supergiants are A(J — 
K) ~0.3 and A{H — K) ~0.1, which each correspond to AA^g ~0.2 using Eqn. ([T]). Hence, if 
the spectral-types of the RSGs are asymmetrically distributed about the mean spectral- type, 
this may produce a derived extinction out by as much as 0.2, consistent with the difference 
between the extinctions derived here and in FMR06. 



3.2.5. Luminosities 

We take the extinctions toward each star derived in Sect. 13.2.41 in conjunction with the 
distance to the cluster estimated in Sect. 13.2.21 to determine the absolute Ks-hand magni- 
tudes of the stars. To convert these to bolometric luminositi es {L-^) we interpolate over the 



contemporary bolometric corrections BCk for RSGs, given in lLevesque et al.l (120051 ). for the 
stellar temperatures derived in Sect. 13.2.31 

The absolute uncertainty in each determination is ^J\5AK + 5BCK + 5Dc\). As we can 
confidently make the approximation that all the stars are all located at the same distance, 
the uncertainty in distance 5-Dci) can be neglected when analysing the luminosity spread of 
the stars to infer the cluster's age (see Sect. 14.11) . The uncertainties in both Ak and BCk 
are carried forward from the error in Tefj, and are determined by substituting the upper and 
lower limits to the stars' temperatures. 

We list the stars' luminosities in Table [2], along with uncertainties which include the 
error in the cluster distance. The newly-derived values are similar to those quoted in FMR06, 
typically within ±0.3dex. As stated above, FMR06 used the same extinction towards all 
stars in the cluster, an approximation which breaks down if there are large variations in 
the interstellar extinction across the field or if a star has extra circumstellar extinction. 
Also, our high S/N, high- resolution spectra give a more accurate picture of variations in CO 
equivalent width, and hence better constrained stellar temperatures - key in evaluating the 
stars' bolometric corrections. For these reasons, we conclude the bolometric luminosities 
derived here to be the more accurate than those quoted in FMR06. 
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3.2.6. Spectral energy distributions 



Using the Galactic plane surveys of MSX and GLIMPSE ( Egan et all200ll : lBenjamin et al 



20031 ). as well as 2MASS, we have collated IR photometry for all stars observed here. For the 
brighter stars, mid-IR photometry is unavailable due to the stars saturating in the images 
(e.g. FOl and F02); while fainter stars in crowded regions (e.g. F15) are dwarfed by brighter 
nearby stars (e.g. F09). We rejected all upper-limit measurements and all detections fainter 
than lOo". 

In Fig. [6]we plot the spectral energy distribution (SED) of each star. Also plotted in the 
figure are the stars' de-reddened photometry, which were calculated using the interstellar ex- 



tinctio n to ward each star in conjuction with the extinction-laws quoted by llndebetouw et al. 



( I2OO5I ) and iMessineo et al.l (120051 ). We have overplotted blackbody curves appropriate for 



the stellar temperatures and luminosities calculated in Sects. 13.2.31 and 13.2.51 

In all cases, the blackbodies provide excellent fits to the near-IR photometry, even in 
the case of F15 where a less accurate method of temperature estimation was possible. This 
serves to validate the stellar luminosities and temperatures of the stars derived above (c.f. 
Fig. 13 of FMR06, where poorer fits to the IR photometry were obtained.). 

In all cases where photometry is available, the mid-IR MSX data shows that the stars 
have considerable excess emission. This is indicative of warm circumstellar dust, a product 
of the high mass-loss rates of the stars. Also, the SEDs appear to show bumps around 12/im, 
which can be understood as silicate emission from the oxygen-rich dust. A detailed study of 
the circumstellar material around these objects will be the subject of a future paper. 



4. Discussion 
4.1. Cluster age 

In Fig. [71 we plot the derived temperatures and luminosities of the cool, luminous stars 
on a H-R diagram. Also pl otted are isochrones taken from the stellar evolutionary models 



of iMeynet &: Maederl (120001 ) which include the effects of rotation, and have initial rotational 
velocity set to 300 km s~^. As we can confidently make the approximation that all the stars 
are at the same distance, we do not include the error in distance on each data-point. The 
magnitude of the error in when this uncertainty is included is shown on the right of the 
plot. 



The figure shows that the temperatures and luminosity spread of the RSGs are well- 
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Fig. 6. — : Spectral energy distributions of cluster members. Crosses are raw photometry from 
2MASS, GLIMPSE &; MSX; filled circles are dereddened according to the stellar extictions derived 
in Sect. I3.2.4[ Over plotted in green is a black-body curve, with T determined from spectral-type 
of star and luminosity according to dereddened K-band magnitude at a distance of 5.88kpc. 
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matched by the 12Myr isochrone. The lOMyr isochrone cannot reproduce the low-luminosity 
stars, while the 14Myr isochrone is too faint to fit the high- luminosity stars. When the 
uncertainty in cluster distance is taken into account (ADci in Fig. [7]), the lOMyr and 14Myr 
isochrones appear more reasonable. We experimented with different evoluti onary models. 



namely non-rotating models with varying ma ss-loss rates and metallicity (jSchaller et al. 



I992I : ISchaerer et al.lll993l : iMeynet et al.lll994l ). We found generally that the non-rotating 
models gave ages that were ;^2Myr younger. We settle on an age estimate for RSGCl of 
12±2Myr. 

We note that the YHG F15 does not lie on the 12Myr isochrone in Fig. [3 Fitting 
the star with the rotating Geneva isochrones, we get an age of lOMyr (pre blue-loop) or 
8Myr (post blue- loop). For the specific models used in this analysis, the masses of RSGs 
in a 12Myr-old cluster do not experience a blue-loop. However, the specific masses of stars 
which experie nce blue- loops are extremely sensitive to the input physics, such as rotation 
(see Fig. 1 of iHirschi et al.l 120041 ). Hence, a blue- loop may be introduced for stellar initial 
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Fig. 7. — : H-R diagra m, showing the locat i on of the RSGs and one YHG. Overplotted are rotating 
Geneva isochrones from iMevnet &: Maeder with ages 10, 12 and 14 Myr. The luminosity 

spread of the RSGs is well-matched by the 12Myr isochrone. Errors on the data-points do not 
include the uncertainty in cluster distance ADd - the magnitude of this error is indicated at the 
right of the plot. The data-point separated from the rest is the YHG F15. 



masses relevant to RSGCl simply by changing the rotational speed. Additionally, it is likely 
that blue-loops would be affected by the inclusion of extra physics (e.g. magnetic fields). In 
summary. We do not necessarily interpret F15's location in the H-R diagram as evidence of 
cluster non-coevality, it could simply be that the input physics of the evolutionary models 
used in the analysis are not fine-tuned to this cluster. 

If the RSGs of this cluster do experience some form of blue-loop, then the position of 
F15 on the H-R diagram is consistent with it evolving away from the RSGs. A post-RSG 
nature for this star would make it a member of a very exclusive club - arguably IRC +10420 
is the only object which is widely acc epted to be a post-RSG , though a case has also been 



argued for HD 179821 (see review of lOudmaijer et al.ll2008l ). We note that F15 does not 



exhibit the same considerable IR excess, nor the bright maser emission of IRC +10420. This 
may be due to F15's lower initial mass - the Geneva models imply ^18Mq for F15, while 
the larger luminosity of IRC +10420 makes it consistent with a star of initial mass ~4OM0 
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(see also Sect. 14.41) . 



The cluster age we derive is slightly greater than the age of ;^9Myr derived in FMR06. 
This previous estimate was determined by comp aring the luminosity spread of the RSGs to 
that predicted by the non-rotating isochones of ISchaller et al.l (119921 ) as a function of age. 
The more rigourous investigation of the stars' luminosities in this present paper results in 
a larger luminosity spread for the RSGs, while the non-rotating models do not reproduce 
the higher luminosities. Much better agreement is found between the contemporary rotating 
models and the new luminosity estimates. 

Finally, we remark that the observed temperatures of the stars are systematically cooler 
than the isochrones (see also Fig. [TT]) . This could be reconciled by increasing either the 
relative metal abundances or the stellar rotational velocities. Both lead to slightly increased 
stellar radii, the former due to the increased opacity of the envelope, the latter due to the 
lower effective gravities. A super-solar metallicity would certainly be consistent with the 
Galactic metallicity gradient and the cluster's Galacto-centric distance (~3kpc). However, 
given the well-known disparities between observations and theory in the field of RSGs, we 
attach a cautionary note to any conclusions derived from this evidence. While recent progres s 
has been made in uniting theory and observation at solar metallicity (Levesque et al.ll2005l ). 
discrepancies still exist at sub-solar metallicities (iLevesque et al.ll2006l ). The location in the 
Galaxy of the Scutum clusters would seem to make non-solar metallicities likely. Accurate 
abundance measurements of the clusters would make them ideal testbeds for evolutionary 
models, as well as probes of the Galactic metallicity gradient. 



4.2. Cluster mass 

To determine the cluster mass we employ the same Monte- Carlo technique used in 
FMR06 and DFK07. We generate a synthetic cluster of a pre-defined initial mass, containing 
stars whose ma sses are random ly drawn from a distribution consistent with a Salpeter initial 
mass function ( Salpetej 1955 ). Then, for a given cluster age, we determine the present-day 



masses, temperatures and luminosities from the Geneva isochrones used in Sect. 14. 1[ We 
then count the number of RSGs in the cluster, where we define a RSG as a star whose 
temperature is lower than 4000K and luminosity greater than IO^Lq. As this is a random 
process, we repeat each simulation 1000 times to reduce statistical noise. 

In Fig. [HI, we plot the number of RSGs contained in a synthetic cluster as a function 
of cluster age, for two initial cluster masses: 2O,OOOM0 and 4O,OOOM0. The plot shows 
that significant numbers of RSGs begin to be seen after ~7Myr - in clusters younger than 
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Fig. 8. — : The number of RSGs as a function of a ge for a coeval cluster of a given initial 
mass, calculated using the rotating Geneva isochrones (iMevnet &: Maedeii l200d ) . The age-range 
for RSGCl, determined from isochrone fitting, is indicated on the plot, as is the number of RSGs 
observed in the cluster. 



this the post-MS stars are massive enough to evolve directly to the WN phase, skipping the 
RSG stage. Peaks in the number of RSGs are reached at ~14Myr and ~20Myr, the dip 
in between is caused by the onset of a blue-loop in the stars' evolution for a narrow range 
of initial masses. The mean luminosity of the RSGs as a function of time decreases, as the 
initial masses of the stars in the RSG zone becomes smaller (see Fig. [TJ and discussion in 
Sect. 14.31) . For a cluster containing 14 RSGs with an inferred age of 12±2Myr (see above), 
we find that the initial mass of the cluster must be somewhere in between these two, implying 
an initial mass of RSGCl of (3 ± 1) x lO^M©. 

Using the velocity dispersion of the RSGs, we can compare this value of the cluster's ini- 
tial mass to the cluster's dynamical mass, under the assumption that the cluster is currently 
in virial equilibrium. The dynamical mass M^y^ is derived using the relation. 
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Fig. 9. — : Cumulative apparent magnitude distribution of the cluster, centred on 18^37"^57.4*, 
-6°52'58.11" (J2000), using the 15 cool luminous stars as tracers of the cluster's spatial distribution. 



where r^p is the half-light radius, is the velocity dispersion, G is the gravitational constant, 
and ?7 is a constant which depends on the stel lar distribution wit h radius, and is typically 
taken to be ~10 (see review in Introduction of iMengel et al.l 120021 ). 

For the velocity dispersion, we find = 3.7 km s~^ after the internal uncertainty in the 
wavelength solution (±lkm s~^) has been subtracted in quadrature. To find the half-light 
radius, we plot the cumulative brightness profile of the cluster using stars FOl-15 as tracers 
of the cluster's spatial distribution (see Fig. [9]). We used a cluster centre of 18'^37'"57.4'*, 
-6°52'58.11" (J2000), the approximate mid-point of the cool stars. We experimented with 
moving the cluster centre by up to 0.5', using different photometric bands, and using the 
luminosities derived in Sect. 13.2.51 rather than the raw photometry. From all these methods 
we found that the cluster half-light radius was stable at 0.8 ± 0.1". At the distance derived 
in Sect. 13.2.21 this gives a cluster size rhp = 1.5 ± 0.3pc. 

Using these values, we find the dynamical mass of RSGCl to be M^y^ = (5 ± 1) ?7/10 x 
W^Mq. Due to the extra uncertainty in the density parameter t], we consider this to be an 
order-of-magnitude estimate only, which compares well to the initial cluster mass derived 
using evolutionary models. 
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Table 3:: Properties of the two Scutum RSG clusters. (1): cluster name; (2): total initial cluster 
mass derived from evolutionary models; (3): cluster mass derived from stellar velocity dispersion; 
(4): cluster age; (5): cluster diameter; (6): cluster distance; (7): distance from cluster to the 
Galactic centre; (8): initial mass range of the RSGs within the clusters. RSGC2 values come from 
DFK07. 



(1) 

Cluster 


(2) (3) 

Minit (evol) Mdyn 

(xlO^Mo) 


(4) 
Age 
(Myr) 


(5) 
d 
(pc) 


(6) (7) 
Dq Dgc 
(kpc) 


(8) 
(Mo) 


RSGCl 
RSGC2 


3±1 5±1 
4±1 6±4 


12±2 
17±3 


1.5±0.3 

o 9+1.2 


6.60±0.89 3.2 

5.83 Q 3.5 


14±2 



4.3. Comparison of the two RSG clusters 

Given that, until recently, the largest number of RSGs in any one cluster was 5, the 
discovery of the two Scutum RSG clusters lying so close to one another is remarkable. After 
applying the same analysis techniques to each cluster, we summarize the physical parameters 
of RSGCl and RSGC2 in Tabled 

The radial velocities of each cluster put them at the tangential point of the Galactic 
rotation curve, with Galacto-centric distances of 3-4kpc. The clusters are separated from 
one another by 0.8^*^0 7 ^pc. Their proximity, combined with their similar ages, suggests that 
they were both formed in a region-wide starburst phase some 10-20Myr ago, and that the 
chemical abundances of their natal material should be similar. 

First-order evidence of uniform metallicity between the clusters comes from the me- 
dian spectral-types of the RSGs, which is M3 in each cluster. The average spectral-type 
of RSGs has been shown to be dependent on environment, shifting gradually to later types 
with increasing metallici ty - averages of K5, Ml, and M2 were found for the SMC, LMC 



and Galaxy respectively (lElias et al.lll985l : iMassey &: Olsenll2003l ). Proposed physical expla- 



nations for this include (a) reduced metallicity leading to lower envelope opacities, increased 
stellar radii and hence to systematically lower effective temperatures; or (b) the effect of 
lower metal abundances on the strengths of the diagnostic TiO lines. Regardless, the aver- 
age spectral-type of M3 for the two Scutum clusters would seem to suggest that they have 
similar, possibly super-solar metallicities. The chemical abundances of the stars in these 
clusters will be the subject of a future paper. 



Both from the analysis of evolutionary models, and the velocity dispersion, we find 
similar masses for each cluster - the factor of two difference in the number of RSGs is caused 
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by the difference in cl uster ages. This may also ex plain the slightly larger size of RSGC2: it 



has been suggested by iBastian fc GoodwinI (120061 ) that a cluster of this initial mass and age 



may be out of virial equilibrium. This arises when the left-over natal material is expelled 
from the cluster by the first SNe explosions of the cluster's most massive stars. This leaves 
the remaining stars with a velocity dispersion larger than that of a virialized cluster, resulting 
in cluster expansion. This may explain (i) the slightly different sizes of the clusters, (ii) the 
apparent lack of any obvious diffuse nebular emission in the GLIMPSE images. Also, we 
note that the dynamical masses of each cluster are slightly higher than the 'evolutionary' 
masses, though the errors on the dynamical masses, in particular the density parameter rj, 
make it difficult to draw any firm conclusions from this. 

The difference in ages between the two clusters implies that the initial masses of the 
RSGs in each must be different; in the younger RSGCl we are seeing stars with larger initial 
masses than in the slightly older RSGC2. To invest igate the initial masses o f the stars in 



each cluster, we again use the evolutionary models of lMeynet fc Maedeii (120001 ). For a given 
isochrone, we find the minimum and maximum initial masses and luminosities of all stars 
with effective temperatures cooler than 4000K, i.e. the temperature range of the RSGs. In 
Fig. [TU] we plot the luminosity and mass ranges of cool stars as a function of time. From the 
luminosity ranges of the stars in each cluster, we then determine the initial masses of the 
RSGs. For RSGCl we find M^,init = ISttMe, and for RSGC2 M^,i„it = 14 ± 2Mq. This 
underlines the potential importance of the RSG clusters to the study of stellar evolution: 
they allow us to study large numbers of RSGs at uniform metallicity as a function of initial 
mass. 



4.4. Masers in RSGs 



One aspect of RSG evolution which the two Scutum RSG clusters allow us to study 
is the onset of the maser-active phase. Masers are often observed in ' extreme' cool stars, 
e.g. the RSGs VY CMa, NML Cyg, S Per (see iRichards fc YateslEJ . and refs therein), 
which are also synonymous with large IR excesses. In the standard picture, these phenomena 
are caused by episodes of high mass-loss, producing large amounts of circumstellar material 
which gives rise to the IR excess. The masers themselves originate in the outfiowing material. 
The SiO 43GHz maser is formed low in the wind at higher temperatures; at larger radii the 
formation of dust grains leads to a depletion of SiO. In the SiO maser-forming region the 
outfiow velocity is low, and henc e SiO masers typ ically have radial velocities similar to 
the stellar systemic velocity, Wgys (jJewell et al.lll99ll ). Higher in the wind the H2O 22GHz 
maser forms, and typically has peaks at many velocities between Vsys ± "^00, where Voo is the 
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Fig. 10. — : The minimum and maximum luminosities (top panel) and initial masses (bottom 
panel) of stars with effective temp eratures cooler than 4000 K in a coeval cluster. Calculated using 
the rotating Geneva isochrones of lMevnet &: Maederl (|200d ). 



wind's terminal velocity (typically around 10-30 km s~^). At larger radii still, H2O is photo- 
disassociated to OH, giving rise to the 1612MHz OH maser. Here the outflow has reached 
its terminal velocity, giving the line-profile its typical double-peaked morphology (where the 
separation of the peaks is 2voo and t he centroid is Wgys)- For a more comprehensive review 



of masers in luminous cool stars, see iHabingl (119961 ). 



It is unclear as to whether the maser stage is one which all RSGs will go through, or 
whether only extreme objects will pass through this phase. Much wo rk has been do ne on the 
detailed physical conditions under which masers form (see review of iHabing] 119961 ). however 
in a simplified picture one may say that the presence of the different masers is determined 
by the wind density, i.e. by mass-loss rate and by pulsations. For an increasing mass-loss 
rate, gradually the critical density will be reached in the formation zones of each transition; 
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Fig. 11. — : H-R diagram of the RSGs in RSGCl, showing those stars which are maser sources. 
A 12M vr model isochone is overplotted, from the rotating Geneva models of iMeynet &: Maedei 

hood ). 



while stellar pulsations will create density-contrasts in the outflow, conducive to population 
inversions. That is to say that we expect to see maser emission from those RSGs which have 
the strongest mass-loss rates and which are pulsationally unstable. 

Mass-loss from RSGs is driven by radiation pressure on dust grains in the outer atmo- 
sphere, and hence depends upon the star's effect ive temperature (cool enough to allow dust to 
form) and luminosity (high radiation pressure) (Ivan Loon et al.ll2005l ). Hence we may expect 
to find masers in high-mass RSGs which have the higher luminosities, while among a coeval 
sample of RSGs one may expect to see masers in those objects furthest along their evolution, 
where their path on the HR diagram takes a sharp upturn (see Fig. [7]). Further, as masers are 
linked to stellar pulsations, we may expect to find masers in unstable stars - i.e. stars with 
high Lj,/M^ ratios that are evo lving closer to the Humphreys-Davidson limit at ~ lO^'^L© 
(IHumphreys &: Davidsonlll979l ). This zone of the H-R diagram is often linked to the so-called 
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'modified' Eddington limit, when the contributions of atomic/molecular transitions to the 
continuum opacity are included when calculating the Eddington luminosity. 

To investigate the presence of masers in RSGs, in Fig. [TT]we plot a HR diagram of the 
RSGs in RSGCl, and indicate those stars which are maser sources. In addition to the OH 
maser ob servations presented h e re, Si O and H2O maser observations of RSGCl were also 



taken by lNakashima fc Deguchil (l2006l ). They found that stars FOl, F02, F04 and F13 were 
spatially coincident with SiO emission, while they concluded that the H2O maser emission 
they found was likely to come from Fl^. 

The figure shows that it is the most luminous stars of the cluster which exhibit maser 
emission. Further, F13, which is the source of SiO, H2O and OH masers, appears to be at 
the point of evolving back toward the blue. The figure therefore appears to support the 
hypothesis that maser emission activates (or becomes strong enough to be observed) in the 
latest stages of RSG evolution, when the star's mass-loss rate is highest and the star becomes 
unstable to pulsations. 

This hypothesis could be tested with a study of the second Scutum RSG cluster, RSGC2. 
This cluster also has many RSGs and so we are again seeing stars in both the earlier and 
later RSG stages. As mentioned in Sect. 14. 3[ the stars in RSGC2 are less luminous and hence 
further from the H-D limit. Though they likely have lower initial masses, they have lower 
L^/M^ ratios. Thus, we may not expect to see maser emission from these stars, especially 
those with lower L^. A comprehensive maser study of the RSGC2 region would be able to 
test at which phase of RSG evolution stars become maser-active, while non-detections would 
place lower-limits to the initial-mass requirements to pass through the maser-active phase. 

Finally, we note that the YHG F15 is not observed to have maser emission, unlike 
the prototype post-RSG IRC +10420. Masers are rarely observed around stars hotter than 
~4000, and the presence of maser emission in IRC +10420's outflow is commonly accepted 
as evidence of its rapid evolution away from the RSG phase. That no maser is observed 
in F15 may be indicative of a lower wind-density while in the RSG phase, due to its lower 
initial mass (~18M0 compared to ~4OM0 for IRC +10420). 



''The large beamsize of the 22GHz H2O maser, 73", overlaps several other stars in the cluster. However, 
as these masers are often observed in stars with maser emission from SiO, and as the central velocity o f 
the H2O maser was consistent with that from the SiO maser emission of F13, Nakashima &: Deguchil ( 20061 ) 
concluded that F13 was likely to be the origin of the H2O maser emission. 
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5. Conclusions 

We have presented a comprehensive investigation into the physical properties of the 
himinous cool stars in the massive cluster RSGCl. Using high-resolution spectroscopy we 
have accurately measured the cluster's radial velocity, derived its kinematic distance, and 
reappraised the stars' temperatures and luminosities. We find a larger luminosity spread 
than in the discovery paper, which is well fitted by a cluster age of 12Myr and cluster mass 
of (3±1) X lO^M©. The mass is similar to that of the nearby cluster RSGC2, and we suggest 
that the difference in the number of RSGs in each is due to the separation in cluster ages, 
with RSGC2 being somewhat older. This implies that the initial masses of the RSGs in each 
cluster are different, which we determine to be ^18Mq for RSGCl and ~14M0 for RSGC2. 
The clusters therefore allow the study of RSG evolution as a function of initial mass, while 
constraining the variable of metallicity. Finally, with new 1612MHz radio observations we 
find compeUing evidence that the OH maser is associated with star F13, and collating recent 
maser observations of the cluster we argue that the maser-active phase is associated with 
stars in the latter stages of RSG evolution. 
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